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Abstract: We present numerical results for the production of a W + W~ pair in association 
with a jet at the LHC in QCD at next-to-leading order (NLO). We include effects of the 
decay of the massive vector bosons into leptons with spin correlations and contributions 
from the third generation of massive quarks. The calculation is performed using a semi- 
numerical method for the virtual corrections, and is implemented in MCFM. In addition 
to its importance per se as a test of the Standard Model, this process is an important 
background to searches for the Higgs boson and to many new physics searches. As an 
example, we study the impact of NLO corrections to W + W~+ jet production on the 
search for a Higgs boson at the LHC. 
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1. Introduction 

The search for the Higgs boson at the Large Hadron Collider (LHC) will rely on analysing 
many types of events related to its different production and decay modes p], E|. If a 
Standard Model (SM) Higgs boson mass lies in the range 155 < ran < 185 GeV, the 
production of a Higgs boson that decays to W pairs, H — > W + W~, is expected to be a 
significant channel. It could even be the discovery mode for the Higgs boson, particularly 
if the mass of the Higgs boson lies very close to the threshold for the production of W 
pairs, Mu = 2M^. Note that the W-bosons can both be real {Mu > 2My/) or virtual 
{M H < 2M W ). 

There are two main mechanisms for producing a Higgs boson that leads to such final 
states at the LHC. The Higgs boson can be produced as the result of gluon fusion, with the 
Higgs boson coupling to an intermediate heavy quark loop (Fig. |l](a)), with further QCD 
radiation leading to additional jets being observed in the final state. Alternatively, the 
Higgs boson may be produced with a sizeable cross section via weak boson fusion (WBF - 
see Fig. |l|(b)), in which the decay products of the Higgs boson are naturally accompanied 
by two forward jets. Therefore the channel W + W~ accompanied by 0, 1 or 2 jets will be 
subject to intense scrutiny at the LHC and all SM backgrounds should be investigated as 
fully as possible. 

If such a Higgs boson exists, the WW + jet sample is expected to contain events 
originating from a Higgs boson that is predominantly produced via gluon fusion, with the 
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Figure 1: The two primary mechanisms for producing a Higgs boson in association with jets at 
the LHC. 



largest background being the normal QCD production of WW pairs. This signal has been 
calculated not only to NLO [|, [| but also up to NNLO f|, |], whilst the NLO corrections 
to this background are also known (7|, ||, [|. 

If instead the sample is selected by requiring two forward jets then the putative Higgs 
boson is preferentially produced via the WBF mechanism. This is now expected to be 
the single most significant channel in the range 130 < mn < 190 GeV (see for example 
Ref. [0]). In this case, the signal has been known to NLO for some time 111], [l2|, [l3]| , whilst 
the NLO calculation of the dominant background from the production of WW + 2 jets is 
currently unknown. 

It has recently been suggested that the Higgs signal significance can be improved 
by considering also the final state which consists of the leptonic decays of W-p&irs, plus 
at least one additional jet [14]. In this study, only a single jet is demanded at large 
rapidity. This has the effect of substantially reducing the backgrounds whilst not unduly 
compromising the signal because of the nature of the WBF events. This renders such a 
search feasible, albeit with significant SM backgrounds whose size should clearly be assessed 
beyond the LO approximation. Dominant among these is the production of Typ^+jet, the 
NLO corrections to which we present in this paper. Since the final state is identified 
through the leptonic decay modes of the W + W~ system to £ + £~ and missing energy, we 
include in our calculation the leptonic decays of the W's and their spin correlations . The 
calculation is performed using a recently-developed semi-numerical method for computing 
virtual 1-loop corrections [16, 17, |lc|. 

The importance of this calculation has already been recognized by its appearance 
on the list of next-to-leading order priorities assembled at the Les Houches workshop in 



In the final stages of preparation of this paper, an independent calculation of this process has been 
presented E3|. The results of that calculation appear consistent with the ones indicated by this paper, 
although Ref. |l5| does not contain an evaluation of the NLO corrections at a specific phase space point 
that we could exactly compare with in a straightforward way. We note however that the calculation of 
Ref. does not include the decay products of the W bosons and their correlations. 
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2005 [19|. This is motivated largely by its obvious importance to the search for the Higgs 
boson and for other signals of new physics at the LHC. Nevertheless, the VFVF+jet process 
is interesting in its own right as a test of our theoretical understanding of the SM. In the 
same way that the SM has been probed experimentally at the Tevatron by measurements 
of the WW and W-\- jets cross sections (for recent examples, see Refs. [2C] and pi]]), the 



rate of events at the LHC will allow similar studies of VFVF+jet events at even higher 
energy scales. In addition, this calculation is a necessary stepping stone to highly-desired 



NLO calculations of 2 — ► 4 processes such as the aforementioned WW + 2 jet process [19], 
as well as W + 3 jet production. 

The structure of the paper is as follows. In section |2] we discuss the overall structure 
of the calculation, with attention paid to the finite contribution arising from internal loops 
of massive quarks. Section || presents the results of our calculation with two types of event 
selection, one indicative of a typical inclusive analysis at the LHC and the other relevant 
for a search for the Higgs boson. To that end, we perform a parton-level analysis of all 
the major background processes using the general purpose NLO program MCFM [B]. Our 
conclusions are contained in section [|. Finally, the appendices consist of analytical results 
for some of the massive triangle diagrams that appear in our calculation, together with 
explicit numerical results for the 1-loop corrections at a particular phase space point. 

2. Structure of the calculation 
2.1 Lowest order 

Let us consider as the lowest order, O («g W a s ), reference process the one in which the up 
quarks annihilate into a W + W~ pair, 

- u(pi) + u(pa) + W+(l(p 3 ) + £(p A )) + W-(£( P5 ) + I(p 6 )) + g(p 7 ) . (2.1) 

The Feynman diagrams and the momentum assignments are shown in Fig. ||. All momenta 
are outgoing. The results for a positive helicity gluon and left handed quark line are 
determined in terms of the two primitive amplitudes given below p2| 2 . 



A (a) _ (15) 
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(15)[56][42](3|2+4|7) <3| 1 + 5|6)(1|2 + 7|4) 
*234 (27) 



(15)(1|2 + 7|6)(3|5 + 6|4) 



(2.2) 



(17)(27)s 56 s 34 si27 
(13)(1|2 + 7|4)(5|3 +4|6) - (35) [46] <1| (3 + 4)(2 + 7)|1) 



(2.3) 



where Sy = [p% +pj) 2 , Sijk = (pi +Pj +Pk) 2 ■ (ij) and [ij] are the standard spinor products 
for massless vectors such that (ij) [ji] = Sj,- . The results for other helicities may be obtained 



2 Note that apart from the overall sign, this agrees with the result obtained in p3] Eqs.(2.22) and (2.23) 
for u{px) + u(p 2 ) + W-(£(p 3 ) + e( Pi )) + W+(£(p 5 ) + t(p e )) + g(p 7 ). 



- 3 - 



Ps 



"Pe 



P7 



Ps 



"Pe 



lP? 



Pv 



Ps <? "Pe 



Pi S -Pa Pi 

P3 > -P4 



"Pa Pi 



Ps ? -P 4 



(a) 



"Pa 
Ps > -P 4 



Ps -P 



Ps -P 4 




Ps -P 



Ps -P 4 




Figure 2: Leading order diagrams for up quark annihilation to W + W + g 



by permutations of the results in Eqs. (2.2, |2.3| ). Since we are interested in W + W _ +jet as 
a background to Higgs boson or new physics searches, we need to go beyond the zero-width 
approximation and in general we have q v ^ My. We introduce the W propagators with a 
Breit-Wigner form, 



Pw{s) 



s - Myy + iT w M w 



(2.4) 



If needed one can also introduce a width for the Z-boson. To do this in a gauge invariant 
(though not unique) way one can multiply the whole amplitude by an additional factor [ 23 



Pz{s) 



s- Mf 



s- Ml + iTzMz 



(2.5) 



To obtain the full amplitude the results of Eqs.(^, |2.3| ) must be dressed with appro- 
priate couplings, 



A A (ui,U2,V3,ei,IJL 5 ,^6,0?) 



A^> (1, 2, 3, 4, 5, 6, 7) + C Lj {u} (s l27 )A\ u > (1, 2, 3, 4, 5, 6, 7) 



C' RjM ( Sl27 )4 6) (2,l,3 J 4,5,6,7) 



(2.6) 



where the couplings Cli u \ and Cri u \ that appear will be defined shortly and H is an 
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overall factor given by, 

H = igs9wPw( s ^)Pw( s 56)Pz(sx27) • (2.7) 
The corresponding results for processes involving d quarks can be written as follows, 
A A (df,df,v 3 ,et,V5,v 6 ,g?) = 

4 a) (1,2,5,6, 3, 4, 7) + C LM} {a m )A® (1, 2, 5, 6, 3, 4, 7) 



Aj(d?,d%,v 3 ,et,^,V6,g?) = V2tf ll2 H 



^, {4 (si 27 )4 b) (2,l,5,6,3,4,7) 



, (2.8) 



where we now specify the left- and right-handed couplings of up- and down-type quarks 

as, 

s(2T 3 / - 2Qsin 2 6 W ) 



2Q sin 2 9 W + 



s-M\ 



C 



R i u \( s ) = ±2(5sin 2 #ty 



(2- 



(2.10) 



In these formulae, Tjj = ±| and is the colour matrix in the fundamental representation 
for the emitted gluon, normalized so that Tr(^i B ) = ^5 AB . g s represents the strong 
coupling and the weak coupling is related to the Fermi constant by g^/{%My^) = Gf/V%- 
The extra sign for the d-type quarks is due to the fact that the triple boson term does not 
require the shift (3 <-» 5, 4 «-> 6) when we go from u to d. If it is included anyway then an 
extra minus sign results. 

The leading order cross-section is obtained by considering all possible crossings of the 
quarks and gluon and by summing the squared amplitudes for all final state helicities and 
averaging over the spin and colour of the initial state partons. At leading order the cross- 
section receives contributions from processes with an incoming quark-antiquark pair or 
with one (anti-)quark and one gluon in the initial state. 

2.2 Real corrections 

The real matrix element corrections to + 1 jet production are obtained by including 

all crossings of results for the two basic processes, 

u(pi) + u( P2 ) + W + (£{ P3 ) + Z{p4j) + W-(£{ P5 ) + l(pe)) + g(p 7 ) + g(p 8 ) , 
- u( Pl ) + u(p 2 ) + W + {i(p 3 ) + £(p 4 )) + W-(£{p B ) + Ifa)) + q'(p 7 ) + q(P8) ■ 

(2.11) 



These tree-level matrix elements have been checked against the results of Madgraph p5| . 
At O (Qg W a 2 ) the cross-section receives contributions from processes with two incoming 
gluons and a qq (qq)-paii as well. Soft and collinear singularities are handled using the 
dipole subtraction scheme [26|. 
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2.3 Virtual corrections 

For the virtual corrections we have to consider one-loop QCD corrections to the tree-level 
processes in figure ||. We have altogether 30 diagrams where the W-bosons are attached 
directly to the incoming quark lines and which result from dressing the diagrams of figure 
|2](a) with virtual gluons. In addition, we have the 6 diagrams of Fig. |3| which have two 
charged bosons attached to a fermion loop. 

There are 11 diagrams which result from dressing the diagrams of Fig. ^b with virtual 
gluons. In addition there are the four diagrams of Fig. |4| which have a Z boson attached to 
a fermion loop and 6 bubble diagrams with a fermion loop which vanish trivially because of 
colour conservation. The full one- loop virtual corrections we consider here involve fermionic 
corrections coming from two massless families circulating in the closed fermion loop and 
from a third massive generation, rrif, ^ m t ^ 0. (We could as well choose to set mj, — ► 
and keep only m< finite). 

The massless contribution to the virtual matrix elements have been computed using the 



semi-numerical approach presented in Refs. [16, 17, 18| based on Davydychev reduction of 



tensor integrals followed by scalar reduction of higher dimensional scalar integrals. Results 
were checked using a second method based on Passarino-Veltman reduction p7|]. The 



method is considerably more efficient than the one used for [g^, 29, ^] since we have now 
a booking of computed integrals not only for higher dimensional scalar integrals, but also 
for tensor ones. 

An additional complication compared to the applications of the semi-numerical method 
considered before is that we have closed fermion traces involving one or two 75 matrices. 
We adopt the t'Hooft-Veltman prescription and split 7^ matrices into a four-dimensional 
part 7^ which anti-commutes with 75 and an e-dimensional part 7^ which commutes with 
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Figure 4: Fermion loop diagrams with Z contributing at one loop. 



75, 

7^ = % + %, {%,75} = 0, [%, 15 ) = 0. (2.12) 

For the massive contributions we have to consider box diagrams with two W bosons 
attached to the loop, shown in Fig. |3| and triangle diagrams where a Z boson is attached to 
the quark loop shown in Fig. ||. We computed the triangle diagrams both fully analytically 
(results are given in App. |A]) and numerically using LoopTools [31]. Box integrals have 



been computed only numerically. Results for these amplitudes, as well as for the massless 
amplitudes, are given in Appendix ^ for one randomly chosen phase space point. 

The whole calculation is incorporated into a private version of the general-purpose 
next-to-leading order code MCFM |8|. The code is flexible enough to allow running with 
on-shell or off-shell bosons. 



3. Results 

In the following we present for results for W + W~ + jet production at the LHC {pp colli- 
sions, = 14 TeV). We use the following electroweak parameters: 

M z = 91.188 GeV , M w = 80.419 GeV , M H = 170 GeV , 
T z = 2.49 GeV , T w = 2.06 GeV , G F = 1.1663910" 5 , 

and furthermore we have, 

m t = 172.5 GeV , F t = 1.48 GeV , m b = 4.62 GeV . (3.1) 

For the calculation of WW+jet at NLO (LO) we use MRST2004f4nlo (MRST2004f41o) 
parton distribution functions corresponding to a s {Mz) = 0.1137 (a s (Mz) = 0.1251) p2| . 
Furthermore we use a four-flavour running and four-flavour evolution of pdfs and the sum 
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Figure 5: Scale dependence for the cross section with cuts I only, as specified in the text. 



over partons in the incoming protons runs over light partons only (d,u,s,c). When we 



consider the Higgs search in section 3.2 we use standard, five-flavour parton densities with 
five-flavour running of the coupling for the signal processes and all backgrounds other than 



VKVK+jet. Specifically, we use MRST2004 |3jj. To define a jet we run the inclusive kr 
algorithm with R = 0.6. In the following we consider two sets of cuts, a more inclusive 
one and one designed to suppress QCD W + W~ + 1 jet events and other backgrounds 
compared to H{— > W + W~) + 1 jet. 

3.1 Inclusive cuts 

We consider first a fairly minimal set of cuts to examine the effect of the NLO corrections 
to the WW+ jet process in an inclusive study. We require at least one jet with, 



P tij i>30GeV, l^il < 4.5 : 



(3.2) 



and then some basic requirements on the decay products of the W's: a minimum missing 
transverse energy and two opposite sign charged leptons with, 

iViss > 30 GeV , P Ml > 20 GeV , P Ma > 10 GeV , | %l(/a) | < 2.5 . (3.3) 

Furthermore we require the leptons to be isolated, i.e. we impose, 



R iih >0A, R jA >0A 



Re lt e 2 > 0-2 • 



(3.4) 



Together, we refer to the cuts in equations Q3.2|), (3.3) and (3^4) collectively as "cuts I". 

In figure |5] we plot the total cross section with cuts I. We fix the renormalization scale 
(or the factorization scale) to be equal to 80 GeV and vary the other scale around this 
value. We see that the NLO correction is sizeable, of the order of 25%, and that this 
is not covered by the scale variation of the LO prediction. While the factorization scale 
dependence is quite mild both at LO and at NLO, the dominant variation is due to the 
renormalization scale. This dependence is reduced by roughly a factor two at NLO. 
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A comparison of the shapes of various kinematical distributions in the T^VF+jet process 
NLO with cuts I. The observable Mt,ww> the transverse mass of the W + W~ pair, is 



Since our calculation is implemented in the parton-level Monte Carlo program MCFM, 
we are also able to examine the effect of the QCD corrections on any infrared-safe observ- 
able. In Figure || we compare the LO and NLO predictions at = = 80 GeV for 
a few examples of distributions that are often important in experimental analyses: the 
rapidity rj ]et of the leading jet, the azimuthal angle Acj) e ^ between the two charged leptons, 
the invariant mass of the two charged lepton system M efl and the transverse mass of the 
W + W _ pair defined as in Ref. Q, 

M T) ww = J(fir + E T ,e„) 2 ~®T + PT,e„) 2 , (3.5) 

where, 

®Tw = yJ&^ + M^ , $t = V^ + M^ . (3.6) 

We have rescaled the LO prediction by the _RT-factor which can be read out of table |l| so as 
to be able to compare the shapes of the distributions. The results indicate that the shapes 
of the LO distributions are mostly unchanged at NLO. 

3.2 Application to the Higgs search 

As mentioned already, the VFVF+jet process is expected to be a significant background to 
the H(—> WW)+}et search channel at the LHC. In this channel there are two mechanisms 
by which the Higgs boson can be produced: 
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(a) Higgs production via weak boson fusion, e.g. 

(3.7) 

U w + w~ 

for which the characteristic signal is two jets at large rapidities and the Higgs decay 
products in the central region, which contains little jet activity. 

(b) Higgs production via gluon-gluon fusion 

9 + a -> 9 + H , 



(3.8) 

w + w~ 



where the basic interaction between the Higgs boson and the gluons is via a top quark 
loop, which is usually treated in the infinite top quark mass limit. Although the QCD 
partons are indicated by gluons, the inclusion of diagrams with gluons replaced by 
quarks is understood. 



The NLO corrections to process (a) are well-known [11, 12, 13] and small, at the level of a 
few percent. In contrast, process (b) receives quite large corrections at NLO |35|]. Both of 
these signal processes are included at NLO in the program MCFM. 

In order to make a more reliable assessment of the viability of this search channel we 
will examine the effect of the inclusion of NLO corrections to both signal and background 
processes when using typical search cuts [H. These cuts are designed to enhance the 
WBF process (a), whilst suppressing the large QCD backgrounds. We will be particularly 
interested in examining the behaviour of the WW+jet process in the search region. The 
background processes which we will consider are as follows. 

(c) QCD production of T^-pairs plus a jet, the process for which we have calculated the 
NLO corrections in this paper, 

q + q^W + W~ +g , (3.9) 
where diagrams related by crossing are understood. 

(d) Associated single top production, 

g + b -> W~ +t , 

I (3-10) 

U w+b 

in which the top quark is produced from a 6-quark in one of the incoming protons. 
The NLO corrections to this process have been known for some time and are typically 



fairly small 37]. They are also implemented in MCFM. 



(e) Top pair production, 

9 + 9 t + t , 

U W~b (3.11) 
■ W + b 
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for which the NLO corrections are known and sizeable |38| . In this study events 
in which two Ws and only one 6-quark are observed are considered part of the 
NLO corrections to process (d). This background assumes that both 6-quarks are 
observed and it contributes because our signal contains one or more jets. Although 
the corrections to the ti process are included in MCFM, this does not include the 
decay of the top quarks. Since we require cuts on the leptons produced in the top 
quark decays, we therefore limit our study of this process to LO only. 

(f) A process analogous to the one which we consider in this paper, Z-pair+jet produc- 
tion, 

q + q -» ZZ + g . 

U vv (3.12) 

• e + e~ 

The NLO corrections to this process are currently unknown and could be computed 
precisely as those for qq — > WW + g presented in this paper. Accordingly, MCFM 
implements currently only the LO contribution. We will see however, that the cross- 
section for this process is very small, so that no NLO corrections are required. 

(g) The production of PF-pairs by EW processes, e.g. 

u + d^ W~W + +d+u . (3.13) 



This process is known at NLO [39], but is not implemented in MCFM. Given that 
its cross-section turns out to be very small, we have calculated only its leading order 
contribution by using the Madgraph package [29] . 



In order to suppress the backgrounds (c)-(g) with respect to the Higgs boson signal 



processes (a) and (b), we follow the strategy of Ref. [14] and apply cuts on four additional 
variables. The shapes of the distributions for both the signal and backgrounds are shown 
in Figure |7], with only the inclusive set of cuts ("cuts I") applied. The signal processes are 
calculated using a putative Higgs mass of 170 GeV, above the threshold for production of 
two real Ws. 

In these figures all the curves have been normalized to unit area to indicate only their 
shape, with the corresponding cross sections given in table |l[ We notice that with cuts 
I the tt process is by far the dominant background. Furthermore, since process (g) has 
only a very small cross section even with this minimal set of cuts, we do not plot the 
corresponding distributions in Figure |?] and do not consider it in the analysis any further. 

From these plots, it is clear that WBF production of the Higgs boson is associated 
with two jets forward in rapidity. Therefore, in addition to our previous selection ("cuts 
I"), we require that the leading jet and the second jet, if present, should be in the forward 
region, 

| % |>1.8 k,|>2.5. (3.14) 

Furthermore, the azimuthal angle between the charged leptons and their invariant mass 
should satisfy, 

(f>i 1 t 2 <1.2 mi 1 £ 2 <75GeY. (3.15) 
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Figure 7: Shapes of the signal and background distributions for the rapidities of the leading and 
sub- leading jets (top- left and top-right respectively), the angle between the two leptons (bottom- 
left) and the invariant mass of the lepton pair (bottom-right). The curves shown are: (a) Higgs 
production by WBF (red); (b) by gluon fusion (black); (c) WW+jet (blue); (d) W + t (green); (e) 
tt (magenta); and (f) ZZ+jet (cyan). The two Higgs production modes (by WBF and by gluon 
fusion) and single top production are calculated at NLO. 



The first of these cuts is a clear discriminator between the two Higgs boson signals and 
all the backgrounds, whilst the second particularly discriminates against the background 
from ZZ+jet and top quark events. We refer to the cuts in equations ( [3.14 ) and ( |3.15| ) as 
"cuts II". 

In figure ^ we examine the scale dependence of the VFH^+jet background cross section 
with the additional cuts applied. There is a very large increase in the cross section going 
from LO to NLO and we see very little improvement in scale dependence. These very 
large corrections are typical of calculations in which new sub-processes enter at NLO. For 
WW+jet, the only contributions at LO are from qq, qg and qg initial states. At NLO, in 
the real radiation term, one has also gg and qq (qq) initial states. In particular, the qq-type 
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Table 1: Cross sections (fb) for signal and background cross sections, using cuts I only, as described 
in the text. 



Process 


ctlo (fb) 


(TNLO (fb) 


(a) H -> WW (WBF) 


22.5 


22.3 


(b) H — > WW (gluon fusion) 


35.9 


64.5 


(c) WIF+jet 


220 


269 


(d) W + 1 


210 


216 


(e) tt 


3090 




(f) ZZ+jet 


36.3 




(g) WW+2 jets (EW) 


9.3 





Cuts l+ll 



12 



^ I I I I I 

\ 




LO |i R -(i, n F -80 GeV 




LO n F =n, |i R =80 GeV 




NLO n R =n, n F =80 GeV - - " 




NLO n, F =n, ^ R =80 GeV - - 


1 1 


i i i 



40 60 80 100 120 140 160 

(i [GeV] 

Figure 8: Scale dependence for the cross section with the cuts I and II. 



processes are very similar kinematically to the WBF Higgs signal. Therefore when cuts 
II are applied, there is a large contribution from the qq initial state since the cuts cannot 
suppress such "signal- like" sub-processes. 

The cross sections for the signal and background processes, after both cuts I and II 
have been applied, are collected in table ^. From the cross sections given in this table we 
see that the cuts have been effective in reducing the backgrounds compared to the signal 
processes. Comparing with table [j], we see that the WBF signal process is reduced by 
about a factor of 2 and the gluon fusion process by a factor of 3.5. In contrast, the huge 
backgrounds from top pair production and single top production have been reduced by 
factors of about 300 and 30 respectively. The ZZ+jet process has been rendered negligible 
as a background by the cut on the invariant mass of the lepton pairs. For the VFVF+jet 
process which we have calculated in this paper, the cuts reduce the cross section by a 
factor of 20 at LO. However, since the NLO correction amounts now to around 70%, the 
overall reduction of the QCD production of WW + 1 jet at NLO is only by a factor of 
about 14. As a result this process is the dominant source of background events in our 
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Table 2: Cross sections (fb) for signal and background cross sections, using cuts I+II, as described 
in the text. 



Process 


olo (fb) 


ctnlo (fb) 


(a) H -> WW (WBF) 


10.6 


10.6 


(b) H — > WW (gluon fusion) 


8.6 


18.0 


(c) WW+jet 


11.7 


20.2 


(d) W + 1 


7.8 


7.6 


(e) tt 


12.7 




(f) ZZ+jet 


0.44 





analysis, although a naive estimate of S/vB from this table is about 1. We note however 
that we have not included NLO effects in the estimate of the top pair background. For the 
totally inclusive top cross section, using our choice of parameters, the enhancement from 
LO to NLO corresponds to a K-factor of approximately 1.5. It is therefore reasonable 
to expect that the top pair cross section would also increase at NLO in the phase space 
region selected by cuts I and II. Therefore we expect the VFTy+jet and tt backgrounds to 
be comparable and a significant source of background events in this channel. 



4. Conclusions 

In this paper we have computed the NLO QCD corrections to the production of W pairs in 
association with a jet at hadron colliders, a calculation highly desired by the experimental 
community [O]. The calculation is performed using a semi-numerical approach and imple- 
mented in the general purpose NLO code MCFM. We find that, for cuts typical of LHC 
analyses, in which a jet is defined with transverse momentum above 30 GeV, the effect of 
the QCD corrections is to increase the cross section by 25% for our default scale choice 
(fJ-F = f-R = 80 GeV). Dependence on the renormalization and factorization scales at NLO 
is mild and is decreased by about a factor of two compared to LO. 

We have also performed a parton level analysis of the impact of these QCD corrections 
on the search for a Higgs boson of mass 170 GeV at the LHC, using the channel H — > 
W + W~+ jet. Including the effect of NLO corrections also in the signal processes, Higgs 
production by gluon fusion and by WBF, we find that the IFH^+jet background is one of the 
dominant backgrounds and is comparable to the WBF signal. In this analysis we have used 
a tighter set of cuts to select the WBF events and to suppress the large QCD backgrounds. 
We find that these cuts enhance the effect of the NLO corrections to WW+jet, which now 
increase the LO cross section by 70% in this region. We conclude that any studies of the 
Higgs boson search in this channel must take into account this significant correction to the 
number of expected background events. 
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Figure 9: Triangle diagrams with closed fermion loops. 
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A. Triangle-Z production diagrams 

In this appendix we give analytical results for triangle diagrams where a Z is emitted from 
a closed fermion loop. Contributions where a Higgs is emitted are included in the signal 
process b), eq. |3^. Results are valid for an arbitrary value of the mass m of the quarks in 



We calculate the triangle shown in Fig. [| where all momenta are outgoing ki+k 2 -\-k 3 = 
and to begin with k\ 7^ 0, k\ 7^ 0. The result for the two triangle diagrams (including 
the minus sign for a fermion loop) is, 



the loop. 




2 cos 9w 



gw 




(A.l) 



where, 



T^P{k x ,k 2 ,m) = [ d n l Tr{y 75 



^7^ z J 



y-m y+fa-m y+tyi + fo-m 




} • (A.2) 



The most general form of T consistent with QCD gauge invariance, 




= k%Y 



2 1 jivp 



= 



(A.3) 



can be written as 



Y^p = F 1 (k x ,k 2 ,m) {^[Y^fafalsK + r Tr[ 7 'W'&75]fc 1 2 } 
+ F 2 (k 1 ,k 2 ,m) {Tr[ 7 V^fe7 5 ]^ + Tr[7V7^i75]A: 2 2 } 
+ F 3 (k u k 2 ,m) (^ + fc£){T¥[7'VM275]} 




(A.4) 
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The results for the coefficients F{ are [[5], 41 , 42 1 , 

Fi(ki, k 2 ,m) = -I w + In + J 2 o 
F 2 (fa, k 2 ,m) = +J i - - ^02 
F 3 (fa, k 2 ,m) = -In , 
F±(fa,k 2 ,m) = 0, 



where, 



dz\dz 2 dz 3 5(1 — z\ — z 2 — z 3 ) 



Az{z\ 



\z x z 3 k\ + z 2 z 3 k\ + z 1 z 2 kj 



(A.5) 



(A.6) 



and k 3 = (fa + k 2 ) 2 . For the particular case at hand we are interested in k\ = e.k 2 = 0, so 
we get a contribution from F\, and in the case of off-shell W's from F3. In this limit we 
obtain, 



Fi(fa,k 2 ,m) 



F 3 (ki,k 2 ,m) 



1 



2fa.k 2 
+ (2 + 7 



2 + Am 2 C (kl, 0, (fa + k 2 ) 2 ,m 2 ,m 2 ,m 2 ) 

\B ((k 1 + k 2 ) 2 ,m 2 ,m 2 ) - B (k 2 ,m 2 ,m 2 



+ 



Fi(fa,k 2 ,m) 

' r 5 ((fci + fc 2 ) 2 , m 2 , m 2 ) - B (kf, m 2 ,m 2 ) 



fa.k 2 



(A.7) 



where, 

C (fc 2 , fcf, (fci + k 2 ) 2 , m\ , ml, m 2 ) 
£o(A;f,mi,m!) 



1 



[q 2 - m 2 ][(q + /ci) 2 - m 2 ,]^ + fa 2 ) 2 - rn 



11 ' 



1 



(Ai 



Z7r^ j [q 2 — m 2 ] [(<? + fci) 2 — m 2 ] ' 

with fa 2 = fa + k 2 . We now make the identification of the momenta fa = p\ + p 2 , k 3 = 
P3 + Pi + P5 +P6, k 2 = pj to calculate the physical process. To get the contribution to our 
process we have to multiply T^ vp given in Eq. (A.l) by the tensor L, dependent on the 
helicity of the (12)-line and the polarization of the gluon £ u (p 7 , h-j) with momentum pj. A 
left handed quark line adds a factor, 



S12 



gs t? li2 (l-W\2-) 



(A.9) 



The decay of a Z splitting to a pair of VF-bosons (which both decay to leptons) adds a 
factor, 

~9w cos °w 
2S 3 4S 56 (S127 - M|) 

x (3 - | 7 a|4-)(5 - | 7 /3|6-) x P w (su)Pw(s 5 6)Pz(si 27 ) 



2g pa ( P 3 + ptf - 2gPP(p 5 + p 6 ) a + g^(p 5 + p 6 - p 3 - p 4 ) p 

(A.10) 



Collecting terms we may define A^p as, 



L tvp(hi,h 7 ) = 8V2g s gl v cos9 w (t A ) lli2 P w (s 34 )Pw(s56)Pz(sm) A^ P (hi,h 7 ) , (A.ll 
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where s^- = fa + Pj) 2 ,s ijk = fa + pj + p k ) 2 and, 

Af,u P (hi,h 7 ) = ——= — — ^-(pi (^1)17^^2(^1)) £u(P7,h 7 ) 

16V2si2S34S 56 (Si27 ~ Mf) 

x <3-|7a|4->(5-| 7 /9|6-> 

x ~2g" a (p 3 +p 4 f - 2gPofa +p 6 ) a + g ^(p 5 + P6 - p 3 _ P4 y] . (A.12) 
So the amplitude with a heavy doublet circulating in the loop is given by, 

Ari(^^3,4,/^6,S? 7 ) = £ r AB X L^ up 

f 

= E 1^2 ^ *4 {2T/}r^A^(/ ll , ^) (A.13) 



<?7 



16?r 2 

Note that with these definitions the overall factor H is the same as in Eq. ( |2.7D except for 
the extra factor of and a sign from 2Tg depending on whether a top or bottom quark 
is circulating in the loop. 

Defining T% = Fi(ki,k2,mt) — Fi(k\, we write out in detail the case h\ = 

-l,h 7 = ±l, 

A tI i{qi,Q2^3,ei,^,P 6 ,gf) = ^^jl^A^-l, +1) 

/ 

= J T l[27] [<1 - |3 + 4 - 5 - 6|7-><35>[46] 

4s 3 4S 56 (si27 - M%) L 

+ 2(5 - |3 + 4|6-)(13)[47] - 2(3 - |5 + 6|4— >(15> [67] 

^3[27] 2 



4S1 2 S34S56(S127 ~ M|) 



(12)( S56 -534)<35)[46], (A.14) 



/ 

-^(17) 



., (7- |3 + 4-5 — 6|2— > (35) [46] 
4s 3 4S 56 (si27 - M|) L 

+ 2(3 — ]5 + 6]4— ) (75) [26] -2(5- |3 + 4|6-)(73)[24] 
^ 3 (17) 2 



4si2S 3 4S 56 (si27 - M| 



[12]( S34 - S56 )(35)[46]. (A.15) 



Switching between the two gluon helicities is given by — flipi, with the flipi symmetry 
defined as, 

fli Pl : 1^2, 3^6, 4^5, (ab) <-> [ab] . (A. 16) 
The cases hi = +1, h 7 = ±1 are simply given by, 

Ari(9f ,92>^> e 4 ^5^e,97 /R ) = Ari(q2,qi,V3,ej,^,v G ,g 7 /R ) . (A.17) 
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Similarly to Eq. ( A. 14 ) we define the contribution from the six massive box diagrams 

as, 

g 2 

•AboxCtfi 1 )^ 1 ) ^e!,/^ ,P 6 ,57 7 ) = Jtf~2 H ^tfin ^box(gi 1 ,92 1 ,^3,e|,/i^,i>6,57 7 ) • 

_ (A.18) 

In the next appendix we give numerical results for A^q^ 1 , q^ 1 , ^3, e\ , uq, gi 17 ) and 
^boxl^i 1 , ^3' e 4 ' ^5 ' ^6) 9? 7 ) f° r a randomly chosen phase space point. 

B. Explicit numerical results 

In this section we present explicit numerical results for one phase space point for the 
independent amplitudes entering the one-loop cross sections. 



We define the one-loop virtual amplitudes analogously to Eqs.(|2Tq, ^q , |A.14| and |A.1S| ), 

e.g. 

167T 2 



At(uL,un,v 3 ,ei,^,v 6 ,g?) = -Ms-^yft^H (B.l) 



4? (1, 2, 3, 4, 5, 6, 7) + C Lt{u} (s m )A$> (1, 2, 3, 4, 5, 6, 7) 

and we have similar definitions for the other amplitudes. 

In table ^| we present results for the colour stripped virtual amplitudes A^j (1, 2, 3, 4, 5, 6, 7) 
and A$(l, 2, 3, 4, 5, 6, 7) for posit ive and negative quark and gluon polarizations at the fol- 
lowing randomly chosen phase space point at the LHC - (p x ,Py,Pz, E) [GeV] (all momenta 
are outgoing) at fiR = 80 GeV: 



Pi 

P2 
P3 
Pi 
P5 
P6 
P7 



0.00000000000000, 0.00000000000000, 1021.22119318758, -1021.22119318758), 
0.00000000000000, 0.00000000000000, -238.714576090637, -238.714576090637), 
-71.5344542606618, -183.877222508616, -3.11006048502754, 197.326337775966), 
-9.92033815503652, -76.1125676676337, 49.0057636944973, 91.0664644166627), 
32.5059044554765, 245.099246845329, -495.644737899924, 553.889863453468), 
64.1786550096635, 124.613643938661, -207.896850811885, 250.736037681104), 
-15.2297670494417, -109.723100607740, -124.860731594604, 166.917065951017). 

The poles of the one-loop unrenormalized virtual amplitudes are given by 

4^(1,2,3,4,5,6,7) = 4 /b (l,2,3,4,5,6,7)^c r /4 e ("^ (B.2) 
+ i f Nc (-s 17 r' + N c (-s 17 r*-±(-s 12 )-*)) . (B.3) 



c 



We renormalize the amplitude by adding the MS counterterm 

A 7)C t(l, 2, 3, 4, 5, 6, 7) = -^& ^U 7 (1, 2, 3, 4, 5, 6, 7) , (B.4) 

where 

llN c -4nfT R , N T(l + e)r 2 (l - e) . , 

b = } — , and c r = 4tt Y V —. ' — ^ - . (B.5 

u 3 , v i r(l-2e) 
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Amplitude 


Cp/e 2 o s /(47r) 


c r/ eQ s/(47r) 


a s /(47r) 


\A^(1 H ,2 L ,3,4, 5,6, 7 L )\ 


- 


- 





| J 4J 7 (1 K ,2 L ,3,4, 5, 6,7 L ) 





0.413542 ■ 10~ 16 


0.116250 ■ 10~ 6 


|AE;(l rt ,2 L ,3,4,5,6,7^) 


- 




0.209219 • 10~ 4 


\A b v7 (l H ,2 L ,3,4:, 5, 6,7 L ) 


0.118557 ■ 10 -3 


0.238602 ■ 10~ 3 


0.473769 ■ 10~ 3 


\A^(1 L ,2 H ,3, 4, 5, 6, 7 L )\ 






0.695412 • 10~ 5 


\A% T (l h , 2 H , 3, 4, 5, 6, 7 L )\ 


0.394067 ■ 10~ 4 


0.793076 ■ 10~ 4 


0.127832 ■ 10~ 3 


\A^(1 L ,2",3,4,5,6,7 L )\ 






0.686719 • 10~ 5 


|j4* 7 (1- £ ',2- k ,3,4,5,6,7 x ')| 


0.389141 ■ 10~ 4 


0.783161 ■ 10~ 4 


0.125011 ■ 10~ 3 


1^(1^,2^,3,4,5,6,7^)1 









|A5 7 (l rt ,2^,3,4,5,6,7 K )| 





0.274681 ■ 10~ 16 


0.216005 ■ 10~ 6 


|A 7 (l- K ,2 J ',3,4,5,6,7 it )| 






0.171644 ■ 10~ 4 


\A b v7 {l u , 2 h , 3,4,5,6,7")! 


0.972651 ■ 10~ 4 


0.195750 ■ 10" 3 


0.352339 ■ 10" 3 


1^(1^,2^,3,4,5,6,7-")] 






0.104711 • 10~ 4 


1^(1^, 2", 3, 4,5,6,7")| 


0.593364 ■ 10~ 4 


0.119417 ■ 10" 3 


0.254739 ■ 10~ 3 


|A 7 (l / ',2 J *,3,4 ) 5,6,7 it )| 






0.102672 • 10~ 4 


|^ 7 (1^,2",3,4,5,6,7")| 


0.581805 ■ 10~ 4 


0.117091 ■ 10 -3 


0.247373 • 10 -3 



Table 3: Numerical results for the Born and virtual amplitudes with massless particles in the loop. 



Amplitude 


a s /(47r) 


|A box (l", 2^,3,4,5,6,7^)] 
\A tli {l R , 2 i ,3, 4,5,6, 7 L )\ 
|A box (l L , 2", 3, 4,5,6,7 L )| 
|A tri (l L ,2",3,4,5,6,7 1 ') 
|A box (l", 2^,3,4,5,6,7")! 
| J 4 tri (l«,2 z ',3,4,5,6,7 i{ )| 
\A hox (l L , 2", 3,4,5,6,7")| 
\A tli (l L ,2 R ,3,4,5,6,7 R )\ 


2.9674010~ 8 
8.4802510" 7 
4.8394610" 8 
1.2007410" 7 
7.3561010- 8 
1.1332510 -6 
1.7482310" 8 
1.1486910" 7 



Table 4: Numerical results for the massive triangle and box contributions to the virtual amplitude. 

In table ||| we give numerical results for the finite contributions to the amplitude coming 
from the massive triangles and boxes. 
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